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ABSTRACT 

Theoretical models of star formation generally assume that bipolar outflows are parallel to the mean 
magnetic-field direction in protostellar cores. Here we present results of A1.3 mm dust polarization 
observations toward 16 nearby, low-mass protostars, mapped with ~2.5" resolution at CARMA. The 
results show that magnetic fields in protostellar cores on scales of ~1000 AU are not tightly aligned 
with outflows from the protostars. If one assumes that outflows emerge along the rotation axes of 
circumstellar disks, then our results imply that these disks are not aligned with the fields in the cores 
from which they formed. 

Subject headings: ISM: magnetic fields - magnetic fields - polarization - stars: formation - stars: 
magnetic field - stars: protostars 



1. INTRODUCTION 

Optical polarization measurements of background stars 
show that magnetic fields are well ordered on large scales 
in the low-densi ty envelopes of mo lecular clouds, as in 
the Pipe Nebula ( Franco et al.|2010 ), which suggests that 
these parsec-scale envelopes are magnetically supported 
("subcritical") . Ultimately, however, ambipolar diffusion 
or turbulence allows the formation of "supercritical" dense 
cores in which gravity overwhelms magnetic support. In 
the simplest axisymmetric case, one then expects the field 
lines to be drawn into an hourglass shape by gravita tional 
collapse, form ing a 1000 AU diameter "pseudodisk" (Galli 
k Shu|[l993| . 

Although the magnetic field may not be strong enough 
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to prevent the formation of a protostar within the pseu- 
dodisk, it can have a significant impact on the accretion 
rate onto the star, and on the formation of a rotationally 
supported circumstellar disk in which planets will form. 
In the limiting case of fiux freezing, the field close to the 
protostar becomes strong enough to brake the rotation of 
the infalling gas completely, p reventing th e formation of 
a rotationally supported disk (Galli et al. 2006J. Even if 



ambipolar diffusion allows the field to escape the central 
core, the redistributed fiux tends to be trapped in a ring 
surroundin g the star, greatly reducing the infall rate (Li 
erar][20TTt 



Magnetic braking is less effective, and disks should form 
more easily, if the rotat ion axis of the clou d is at an angle 



to the magnetic field (Joos et al. 2012 1. Indeed, from 
an observational perspective it is clear that circumstellar 
disks typically are not well-aligned with the parsec-scale 
magnetic fields in the surrounding molecular cloud. If 
they were, the bipolar outflows and jets that emerge 
along the axes of these disks would all be p arallel with 
the ambient B-fi eld, but this is not observed: |Menard"^ 



Duchene (2004) have shown that the optical jets from 
classical 1 i'auri stars in the Taurus- Auriga molecular 
cloud are randomly oriented with respect to the parsec- 
scale magnetic field in this cloud; ,Targon et al.^ (20lf| 
obtained a similar result for 28 regions spread over the 
Galaxy, although they do find some evidence for alignment 
of jets from younger. Class and Class I, protostars. 

Polarization observations of background stars are un- 
able to probe the magnetic field morphologies inside the 
dense cores where circumstellar disks form; even at in- 
frared wavelengths the extinction through these regions 
is too great. Mapping the polarized thermal emission 
from dust grains at millimeter and submillimeter wave- 
lengths is the usual means of studying the magnetic fields 
in these regions. Under most circumstances, spinning 
dust grains are expected to align themselves with their 



long axes perp endicular to the magnetic field (Hoang & 



Lazarian 2009 ) , so normally the thermal radiation from 



these grains is weakly polarized perpendicular to the field. 
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Dust polarization maps of many sources have been 
obtained at submillimeter wavelengths with single-dish 
instruments, e.g. with the SCUBA polarimeter on the 
JCMT at 850 fim. ([Matthews et al.|[2009), and with the 



Hertz polarimeter on the (JSU at 350 fim (Dotson et al. 
2010). These maps have angular resolutions ot about 
20" , corresponding to scales of 3000-8000 AU in nearby 
molecular clouds. 

Higher angular resolution is required to study the field 
geometry in the innermost regions of the cores where 
circumstellar disks form. Thus far interferometric polar- 
ization maps with resolutions of a few arcseconds have 
been published for about a dozen sources. Most of this 
work has focused on detailed analyses of in dividual ob- 
jects: examples incl ude NGC 1333- IRAS 4 A ( Girart et al. 
2006 )) IRA S 16293 (| Rao et al.|2009| , and Orion KL ( |'l|aSg 
et ai |[2010 |. In this limited set ot sources, outflows often 
were round to be skewed with respect to the inferred 
magnetic field directions, hinting that circumstellar disks 
may not be tightly aligned with the magnetic fields on 
-1000 AU scales. 

The 1 mm dual-polarization receiver system at CARMA 
(the Combined Array for Research in Millimeter-wave As- 
tronomy) allows us to map the dust polarization toward 
many more sources. Here we present results from the 
TADPOlj^ survey a CARMA key project to study dust 
polarization in protostellar cores. This paper focuses on 
results from nearby, low-mass protostars. We compare 
the field directions inferred from dust-polarization mea- 
surements, with outfiow directions, which indicate the 
axis of the rotationally supported disk. On the '^lOOO 
AU scales probed by our data, magnetic fields appear to 
be randomly oriented with respect to outfiows. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations were made with CARMA between May 
2011 and October 2012. We selected sourc es from catalog s 
of young stellar objects , includ ing J0rgen sen et al. ( 2007 ), 
Matthews et al.| ( |2009[ ), [Tobin et al.| |2010| , and [Enoch 
et al. 1 21)1 1| . We tocus on Class and Class I objects at 
distances of d < 400 pc that are known to have bipolar 
outflows. 

The polarization system consists of dual-polarization 
receivers that are sensitive to right- (R) and left-circular 
(L) polarization, and a spectral-line correlator that mea- 
sures all four cross polarizations (RR, LL, LR, RL) on the 
105 baselines connecting the 15 antennas. The receivers 
comprise a single feed horn, a w aveguide circular polarizer 
(Plambeck & Engargiola"2010|, an orthomode transducer 
(!^lavarrini & Flambcck 2006|, and two mixers. The re- 
ceivers are double-sideband; a phase-switching pattern 
applied to the local oscillator (LO) allows signals in the 
lower (LSB) and upper sidebands (USB) to be separated 
in the correlator. 

For these observations, the correlator was set up with 
three 500 MHz wide bands to measure the dust contin- 
uum, and one 31 MHz wide band to map bipolar outflows. 
The frequency of the flrst LO was 223.821 GHz. The 
three 500 MHz continuum bands were centered at 6.0, 
7.5, and 8.0 GHz in the intermediate frequency (IF). The 
corresponding RF frequencies are equal to the difference 



(LSB) or the sum (USB) of the LO and the IF. The 
spectral-line band was centered at IF = 6.717 GHz, allow- 
ing simultaneous observations of the SiO(5-4) line in the 
LSB, and the C0(2-l) line in the USB. For the spectral 
line measurements, the channel spacing is 0.21 km s^^. 

In addition to the usual gain, passband, and flux cali- 
brations, polarization observations require two additional 
calibrations: "XYphase" and leakage. The XYphase cali- 
bration corrects for phase differences between the L and 
R receivers, and is done by observing an artiflcially po- 
larized noise source of known position angle. The leakage 
corrections compensate for cross-coupling between the L 
and R receivers, and are calibrated by observing a strong 
source (usually the gain calibrator) over a range of paral- 
lactic angles. There are no moving parts in the CARMA 
dual polarization receivers, so the measured leakages are 
stable with time. A typical antenna has a band-averaged 
leakage amplitude (i.e., a voltage coupling from L into R, 
or vice versa) of 6%. 

We perform calibration and imaging with the MIRIAD 
data reduction package (Sault et al. 1995). Using multi- 
frequency synthesis and natural weighting, we create dust- 
continuum maps of all four Stokes parameters (/, Q, U, V). 
The typical beam size is 2.5", which corresponds to a 
resolution of 750 AU at a distance of 300 pc. We produce 
polarization position-angle and intensity maps from the 
Stokes /, Q, and U data, where the position angle of the 
incoming radiation is x = 0.5 arctan ([/ /Q), and the bias- 
corrected polarized intensity is Pc = \/ + — 
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( Vaillancourt 2006). up is the median rms noise in the 
Stokes Q and <7 maps . 

In good weather cp « 0.4 mJy/beam for a single 6- 
hour observation, and can be as low as —0.2 mJy/beam 
when multiple observations are combined. We consider 
it a detection only if Pc > Sep , and the location of the 
polarized emission coincides with a detection of total 
intensity / > 2ai, where ai is the rms noise in the Stokes 
/ map. 

We also generate maps of the red- and blueshiftcd CO 
and SiO line wings to measure outflow directions. We 
generally use CO maps to measure the outflow direction, 
but occasionally we use SiO if the CO emission is too 
complex. We do not attempt to measure polarization in 
the spectral line data because of flne scale structure in 
the polarization leakages. 

3. RESULTS & ANALYSIS 

Of the 27 TADPOL sources within —400 pc, we detected 
dust polarization toward 16, which we focus on in this 
paper. The full set of TADPOL results will be presented 
in a separate paper. 

Figure [l] shows some example results. In all of the plots, 
the dust-polarization vectors have been rotated by 90° 
to show the inferred magnetic fleld direction, and are 
plotted at the Nyquist spatial frequency (two vectors per 
synthesized beam). 

Table [T] lists the results for the 16 TADPOL sources, as 
we ll as tor IR AS 16293, which was previously published 
by Rao et al. (2009). Two of the sources each have two 
distinct outflows, which we consider as independent data, 
thus making a total of 19 entries. Note that NGC 1333- 
IRAS 4A, one of the sources in t he TADPOL sampl e, has 
been mapped in detail before (Girart et al. 2006), and 
was included in our survey as a cross-check. 
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Figure 1. Sample maps from two sources in the TADPOL survey. Line segments show the inferred magnetic field direction; they have 
been rotated by 90° relative to the polarization direction. Vector lengths are proportional to polarized intensity, not fractional polarization. 
The scale bar shows the peak intensity of polarized emission in mjy bm~^. The solid ellipses show the synthesized beam width (FWHM) 
and position angle. Top: Red- and blue-shifted line wings in CO(2-l) for IRAS 4B, and in SiO(5-4) for Ser-emb 8, show the bipolar outflows. 
Velocity ranges are labeled in the figures. Contour levels are in steps of 6cr, beginning at it4(To, where (Tq is the rms noise measured in 
each outfiow map. Bottom: Dust continuum (Stokes I) contours and inferred B-field vectors toward the same two sources. Contours are 
—2, 2, 3, 4, 6, 8, 11, 16, 23, 32, 45 x ctj, where ctj is the rms noise measured in the dust continuum maps, cj = 26 mJy bm~^ for IRAS 4B; 
o"! = 3.3 mJy bm~^ for Ser-emb 8. 



To estimate the outflow direction, we measure the posi- 
tion angles of lines connecting the center of the continuum 
source and the intensity peaks of the red and blue outflow 
lobes; we take the average of the two position angles as 
the outflow direction. As a crude estimate of the un- 
certainty CTo, we use half the difference of the position 
angles. The B-field direction is calculated by averaging 
the B-fleld position angles, weighted by the Stokes I total 
intensity. This gives the fleld direction in the densest 
part of the core. The uncertainty in the B-fleld ctb is the 
average of the B-fleld position angle uncertainties, also 
weighted by the Stokes I total intensity. The total uncer- 
tainty in the angle between the B-fleld and the outflow is 



Figure |2] shows the cumulative distribution function 
(CDF) of the angles between the B-fleld and outflow axes 
from Table [T] The B-fleld and outflow position angles 
we observe are projected onto the plane of the sky. To 
determine if the large scatter in position angle differences 
could be due to projection effects, we compare the results 
with a Monte Carlo simulation. The simulation randomly 
selects pairs of vectors in three dimensions that are within 
20° of one another, and then projects the vectors onto the 
plane of the sky and measures their angular differences. 
The resulting CDF is shown by the dotted curve in Figure 

t Projection effects are not as problematic as one might 
ink: to have a projected separation larger than 20°, the 
two vectors must point almost along the line of sight. 
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A Kolmogorov-Smirnov (K-S) test shows that the proba- 
bility that our data were drawn from the same population 
as the model, where outflows are within 20° of the field, 
is 3 X 10~^, suggesting that the outflows and B-fields in 
our sample are not tightly aligned. 

By comparison, the solid line in Figure |2] shows the 
CDF expected if outflows and B-fields are randomly ori- 
ented. In this case, the CDF is a straight line: all angular 
differences are equally likely. A K-S test shows that the 
probability that the results were drawn from this popula- 
tion is 0.64. 

4. DISCUSSION 

Most analytical models of star formation assume that 
the rotation axis of a protostellar core, its mag netic field 
axis, and its outflow directio n all are parallel ( Shu et al. 
20001 IKonigl & Pudritz"2000 ). The TADPOL results show 



that this is not the case, and find instead that magnetic 
fields on scales of ~1000 AU are not tightly aligned with 
protostellar outflows. 

We consider a number of possible explanations for the 
misalignment of B-fields and outflows. First, since we 
detect a net polarization in the sources in our sample, 
the dust grains are indeed aligned. But are they always 
aligned perpendicular to the magnetic field? If grains are 
mechanically aligned by outflows, then their polari z a.tion 
is expected to be parallel to the B-field ( Gold| 1952 1. Rao 
et al. ( 1998 1 explored this possibility to explain a 90 flip 
in the polarization toward the BN-KL region in Orion. 
If grains are aligned by r adiative torques (e.g. [Lazarian] 
2007 Hoang fc Lazarian||2009 1 , then under some condi- 
tions they can also be stably aligned parallel to the B-field. 
There is possible evidence for anomalous grain alignment 
in at least one of our sources: in Ser-emb 8 (see Figure 
fTI), our results match the larger-scale polar ization from 
^CUBA observations ( [Matthews et al.|2009]) in the north- 
ern core, but in the southern core the B-tields inferred 
from CARMA and SCUBA are perpendicular. While in 
principle one would expect to see a bimodal distribution 
of angles if grains were aligned either perpendicular or 
parallel to the B-field, in fact our observations integrate 
over regions that could contain grains of various orienta- 
tions. This averaging, which occurs both in the plane of 
the sky and along the line of sight, would tend to lead to 
a random distribution of orientations. 

Second, the B-field could be wrapped up azimuthally 
by core rotation. Fo r example, numerical simulations by 
Machida et all ( 2006 1 show that in sources with both weak 
magnetic fields and strong rotation, the field eventually 
becomes perpendicular to the outflow due to azimuthal 
wrapping. As in the case of grain alignment, we might 
expect that a sample of sources with either wrapped 
or straight fields would lead to a bimodal CDF. Again, 
however, averaging over the synthesized beam and along 
the line of sight could lead to a random distri bution. 
Fina lly, s ome simulations, such a s thos e by [Joos et al] 



(20121 and Hennebelle & Ciardi p009|, suggest that 



misalignment of the B-held and the core ro tation axis can 



aid in the formation of circumstellar disks. Li et al. ( 2011 1 
find that disk formation is suppressed in the case where 
fields and core rotation axes are parallel, even when non- 
ideal MHD effects are considered. Hence, these models 
suggest that misalignment may be a necessary condition 
for the formation of disks. 



5. SUMMARY 

Results from the TADPOL survey show that magnetic 
fields on scales of ^1000 AU are not tightly aligned with 
protostellar outflows. If one assumes that outflows emerge 
along the rotation axes of circumstellar disks, then our 
results imply that these disks are not aligned with the 
fields in the cores from which they formed. 

It could be fruitful to investigate whether alignment cor- 
relates with environmental factors, multiplicity, outflow 
velocity, or age. In addition, higher resolution observa- 
tions with ALMA will test these correlations at the ^100 
AU scale of circumstellar disks. 
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Projected angle between outflow and B-field (deg) 

Figure 2. Cumulative distribution function (CDF) of the (projected) angle between the mean magnetic field and outflow axes of the 
sources m Tabled The dashed curve is the CDF from a Monte-Carlo simulation where outflow and B-field directions are oriented within 
20° of one another. The straight line is the CDF for random orientation. 



Table 1 

Observations 



Source 
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Xb (o-b) 








-B (o-Q- 


b) 




d 


Distance 




(J2000) 


(J2000) 


n 




( 






(°) 




(") 


(pc) 


ref.d 


L1448 IRS 2 


03:25:22.4 


30:45:13.2 


139 


(9) 


134 


(5) 


5 


(10) 




3.63 


230 


1 


L1448N(B) 


03:25:36.3 


30:45:14.8 


31 


(6) 
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(2) 


66 


(6) 




2.04 


230 


1 


NGC 1333-IRAS 2A'^ 


03:28:55.6 


31:14:37.1 
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(6) 


16 
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1 
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21 
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NGC 1333-IRAS 4A 


03:29:10.5 


31:13:31.3 
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40 


(10) 
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2 
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54 
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172 
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Note. — Coordinates 


are fitted positions 


of dust emission 


peaks 


The 


outflow angle 


Xo aiid infcrr 


2d magnetic-field 


angle xb 


are measured 



counter-clockwise from north. The angle difference &o-b between the outflow and the B-field is always between 0-90°. The B-field direction is 
assumed to be perpendicular to the direction of the polarized dust emission, d is the distance to the source, ^bm is the geometric mean of the major 
and minor axes of the synthesized beam. 

^ Source has two outflows. 

^ Coordinates from Takahas hi et al.| | |2009| l 
^ Results from^ Rao et al.j (j2009[|. 

^ Distance reference s. 1: |Hirota et al.|i201l [. 2: |Looney et al.|||2000[ . 3: |j0rgensen et al.|||2007| |. 4: |Menten et al.|]2007| . 5: |Loinard et al.|]2008] . 
6: |Imai et"ar] | |2007^ . 7: |Dzib et al.| | |2010| ) 8: [Looney et al.|]2007[ . 9: [Launhardt et al.| | |20'To| k 
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